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ABSTRACT 

Sun-like stars are thought to be regularly disrupted by supermassive black holes (SMBHs) within galactic 
nuclei. Yet, as stars evolve off the main sequence their vulnerability to tidal disruption increases drastically 
as they develop a bifurcated structure consisting of a dense core and a tenuous envelope. Here we present the 
first hydrodynamic simulations of the tidal disruption of giant stars and show that the core has a substantial 
influence on the star's ability to survive the encounter. Stars with more massive cores retain large fractions 
of their envelope mass, even in deep encounters. Accretion flares resulting from the disruption of giant stars 
should last for tens to hundreds of years. Their characteristic signature in transient searches would not be 
the r^^^ decay typically associated with tidal disruption events, but a correlated rise over many orders of 
magnitude in brightness on months to years timescales. We calculate the relative disruption rates of stars of 
varying evolutionary stages in typical galactic centers, then use our results to produce Monte Carlo realizations 
of the expected flaring event populations. We find that the demographics of tidal disruption flares are strongly 
dependent on both stellar and black hole mass, especially near the limiting SMBH mass scale of ~ lO^M©. 
At this black hole mass, we predict a sharp transition in the SMBH flaring diet beyond which all observable 
disruptions arise from evolved stars, accompanied by a dramatic cutoff in the overall tidal disruption flaring rate. 
Black holes less massive than this limiting mass scale will show observable flares from both main sequence and 
evolved stars, with giants contributing up to 10% of the event rate. The relative fractions of stars disrupted at 
different evolutionary states can constrain the properties and distributions of stars in galactic nuclei other than 
our own. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: nuclei — hydrodynamics — 
methods: numerical — stars: evolution — stellar dynamics 



1. INTRODUCTION 

Quasars are rapidly growing black holes lit up by the gas 
they accrete. They are the most dramatic manifestation of the 
more general phenomenon of active galactic nuclei (AGN) 
and they are among the most energetic objects in the uni- 
verse. Several billion years after the Big Bang, the universe 
went through a quasar era when highly-luminous AGN were a 
standard feature of most massive galaxies (Kormendy & Rich- 
stone 1995; Richstone et al. 1998). Since that time, AGN have 
been dying out and the only activity that still occurs in many 
nearby galactic nuclei is weak (Schawinski et al. 2010). All 
that should remain in the centers of local galaxies are the rem- 
nants of quasar-era exponential growth: quiescent SMBHs, 
now dim and starved of fuel (Ho 2008). It is, therefore, not 
surprising that definitive conclusions about the presence of 
local SMBHs are typically drawn from very nearby galaxies 
with little to no AGN activity (Gebhardt et al. 2000; Ferrarese 
& Merritt 2000). The centers of these galaxies are well re- 
solved, revealing the region where the black hole dominates 
the stellar dynamics. 

The question then arises of whether the postulated presence 
of SMBHs lurking in the centers of most galaxies is consistent 
with the apparent quiescence of their nuclei. Quiescent black 
holes are, in fact, nearly black; they may only be lit up by the 
luminance of accreting matter. We do not directly know how 
much gas there is near most SMBHs, and there is no a priori 
reason why nuclear regions should be swept completely clean 
of gas. The distribution of stars in dense clusters that sur- 
round SMBHs, on the other hand, is much better constrained. 
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These densely packed stars trace complicated and wandering 
orbits under the combined influence of all the other stars and 
the black hole itself. If a star wanders too close to the black 
hole it is violently ripped apart by the hole's tidal field (eg, 
Hifls 1975; Frank 1978; Rees 1988). About half of the de- 
bris of tidal disruption eventually falls back and accretes onto 
the SMBH. This accretion powers a flare which is a defini- 
tive sign of the presence of an otherwise quiescent SMBH 
and a powerful diagnostic of its properties (Rees 1988). Tidal 
disruption events are expected to be relatively rare, on the 
order of one per ~ 10^ years per galaxy (e.g. Magorrian & 
Tremaine 1999; Wang & Merritt 2004). Depending primar- 
ily on the structural properties of the disrupted star, an ultra- 
luminous transient signal could persist steadily for months to 
at most tens of years; thereafter the flare would rapidly fade. 
In a given galaxy, this luminous flaring activity would then 
have a short duty cycle. Quiescent SMBHs should greatly 
outnumber active ones. Observational constraints would not, 
therefore, be stringent until we had observed enough candi- 
dates to constitute a proper ensemble average. However, the 
long decay tails of these flares may account for an appreciable 
fraction of the total low-luminosity AGN activity in the local 
universe (Milosavljevic et al. 2006). 

The critical pericenter distance for tidal disruption is the 
tidal radius. 



where and are the stellar mass and radius, and Mbh is 
the black hole mass (Hills 1975). The tidal radius is larger 
than the black hole's horizon, = 2GMbh/c^, for solar type 
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stars as long as the black hole is less massive than about 10^ 
solar masses. In encounters with more massive black holes, 
solar type stars may pass the horizon undisrupted and are ef- 
fectively swallowed whole. Such events would leave little 
electromagnetic signature (although a portion of the gravita- 
tional wave signature would remain, Kobayashi et al. 2004). 
However, for a given stellar and black hole mass, evolved stars 
have larger tidal radii than main sequence (MS) stars and are 
therefore more vulnerable to tidal disruption. Furthermore, 
giant branch stars are the only stars that can produce observ- 
able tidal disruption flares in encounters with the most mas- 
sive black holes > lO^M©. 

Motivated by these facts, we examine the importance of 
stellar evolution in the context of tidal disruption. A sun-like 
star spends the majority of its lifetime on the MS, Tms ~ 10^^ 
years, followed by a relatively brief period of post-MS (giant 
branch) evolution, Tg ~ 10^ years, once its central supply of 
hydrogen fuel is exhausted. As nuclear reactions slow, the 
stellar core loses pressure support, and, in approximately a 
thermal diffusion time, the star ascends the giant branch as its 
outer layers expand in response to the core's collapse. Giant- 
branch stars are much less common than MS stars in a typical 
stellar population due to the ratio of their lifetimes, which for 
a solar mass star is Tg/rms ~ 10~^. However, their large radii 
imply that they are exceptionally vulnerable to tidal disruption 
during these brief periods. The contribution of giant stars to 
the tidal disruption event rate and the luminosity function of 
AGN will depend on the competing effects of their enhanced 
cross- section and their relative rarity. 

Understanding the details of giant star disruption events and 
their contribution to the SMBH flaring population is the focus 
of this work. We use several methods to study this problem. In 
Section 2, we discuss the calculation of detailed stellar evolu- 
tion models and outline the importance of stellar evolution in 
the context of tidal disruption. The non-linear dynamics of the 
encounters themselves must be understood through hydrody- 
namic simulations; this is particularly true for post-MS stars 
which are not well described by a simple single-polytrope 
model. In Section 3, we describe how we derive giant star ini- 
tial models from our stellar evolution calculations, our meth- 
ods of hydrodynamic simulation, and the results of our sim- 
ulations of close encounters between giant-branch stars and 
SMBHs. In Section 4, we calculate the rates of tidal disrup- 
tion that result from the two-body relaxation driven random 
walk of nuclear cluster stars in angular momentum space. We 
focus on the relative rates of disruption of stars in different 
evolutionary states. In Section 5, we combine our stellar evo- 
lution, hydrodynamic, and rate calculations of Sections 2-4 
and present Monte Carlo realizations of flaring events. We 
discuss the demographics of tidal disruption-powered flaring 
events as a function of black hole mass, the contribution of 
giant stars to the luminosity function of local AGN, and the 
detection of flares due to the disruption of giant stars. 

2. STELLAR EVOLUTION IN THE CONTEXT OF TIDAL 
DISRUPTION 

Stellar evolution naturally enriches the physics of tidal dis- 
ruption and leads to a large diversity of accretion flare events. 
The effects of stellar evolution are most pronounced as stars 
evolve off the MS. Their radius can change by orders of mag- 
nitude, they often suffer some degree of mass loss, and the 
density contrast between the core and the envelope dramati- 
cally increases. In this section, we will outline the basic char- 
acteristics of post-MS stellar evolution and provide some in- 



tuition for how the capture rates and properties of the tidal 
disruption events might change as a stars evolve off the MS. 

2.1. Stellar evolution models 

To capture the intricacies of post-MS stellar evolution we 
use the open source Modules for Experiments in Stellar As- 
trophysics (MESA) stellar evolution code (version 3290, Pax- 
ton et al. 2011). MESA solves the stellar evolution equations 
within the Lagrangian formalism, closed by a tabulated equa- 
tion of state which blends OPAL, SCVH, PC, and Helmholz 
equations of state to treat a wide range of fluid conditions 
and ionization states (Paxton et al. 2011; Rogers & Nay- 
fonov 2002; Saumon et al. 1995; Potekhin & Chabrier 2010; 
Timmes & Swesty 2000). MESA handles short evolutionary 
periods (for example, a helium core flash) by switching to 
a hydrodynamic solution which allows non-zero velocities at 
cell interfaces. This capability gives MESA the ability to ef- 
ficiently evolve stars from their pre-MS collapse all the way 
to the formation of a white dwarf. For our purposes this is 
essential because it allows us to capture the entire lifetime of 
stars. 

We compute models of solar metallically stars in the mass 
range of 0.95 -SM©. To do so, we employ an extension to 
mesa's basic hydrogen and helium nuclear burning network 
called agb . net. Stellar mass loss is modeled for red giant 
(RG) stars using Reimers's formula, 

«» = 4.>0-.V(^)(|l)(^)"'«oyr-'. (2, 

with % = 0.5 (Reimers 1975). Mass loss on the asymptotic 
giant branch (AGB) is given by the Blocker formula, 

«, = 1.932xl0-,.(^^)"(|^)(^)"''«e,r-', 

(3) 

with rjB =0.1 (Bloecker 1995). A switch between the wind 
schemes is triggered when the central hydrogen and helium 
are depleted. The efficiencies of wind mass loss rates are rel- 
atively uncertain (e.g. Habing & Olofsson 2003), and they af- 
fect some properties of giant stars. The efficiencies % = 0.5 
and T^B = 0. 1 are typical values used in the fiducial IMq mod- 
els presented with the MESA code (Paxton et al. 201 1). 

2.2. Post-MS stellar evolution 

Stars exhaust their central reserves of hydrogen fuel in 
about 10^^(M^/Mq)-^-^ years (Hansen et al. 2004). Af- 
ter this, they begin to evolve off the MS. Nuclear reactions 
slow in the stellar core, which contracts under the influence of 
its self-gravity as it loses pressure support in approximately a 
thermal diffusion timescale. Hydrogen burns to helium off- 
center in a thick shell that continues to add mass to the now 
inert helium core (Iliadis 2007). The star's structure reacts 
to the core's contraction through an expansion of the star's 
outer layers. During this phase, the star's luminosity increases 
greatly, while its effective temperature decreases causing it to 
ascend the RG branch in the Hertzsprung-Russell (HR) dia- 
gram. The star develops an increasingly bifurcated structure 
with a dense, radiative core and a convective envelope. In 
stars less than about 2.2M0, the helium core collapses to the 
point where it becomes degenerate. When helium burning 
eventually ignites it does so in degenerate material, causing a 
dynamic event known as the helium core flash. Above about 
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Fig. 1. — Radius (left) and mean density (right) evolution of post-MS stars. In the left panel, stages are marked as: main sequence (MS), red giant (RG), 
horizontal branch (HB), and asymptotic giant branch (AGB). In the right panel, we plot the mean density of stars between 0.95 and 5 Mq as a function of 
fractional lifetime. The mean density describes the stars' vulnerability to tidal disruption (1). Periods of low density occupy an increasing portion of the lifetime 
of more massive stars. 



2.2M0, helium ignition occurs off-center in non-degenerate 
material, and the star does not reach as extreme radii while on 
the RG branch (Maeder 2009). 

Following helium ignition, the core expands and the star 
stably burns helium to carbon and oxygen for a period of 
about Thb ~ 10^ years along the horizontal branch (HB). Dur- 
ing this time the core mass grows significantly as hydrogen 
continues to burn to helium in a shell surrounding the core 
(Maeder 2009). Finally, the star may ascend to the tip of the 
giant branch one final time as its carbon oxygen core becomes 
degenerate. Stars above (about) 9Mq will go through addi- 
tional burning phases (Maeder 2009). For the intermediate 
mass stars considered here, at the tip of the AGB, the stars' 
luminosity becomes so great that the bulk of the envelope 
is driven off via a combination of strong winds and thermal 
pulses, eventually exposing the bare carbon-oxygen core - a 
proto- white dwarf (Habing & Olofsson 2003). 

In the left panel of Figure 1, we show typical evolution- 
ary stages for stars with \AMq and 2.^7 Mq zero age MS 
(ZAMS) mass. For simplicity we define here the transition 
from the MS to the sub-giant portion of the RG branch as 
being when a star first exceeds 2.5 times its ZAMS radius. 
The post-MS evolution is a very sensitive function of the ini- 
tial stellar mass. In particular, the peak radii and timescales 
of the RG and AGB phases vary considerably with mass. 
A common feature is that while on the giant branch, stars 
with a wide range of masses spend considerable time hav- 
ing radii in the range ~ 5 - 20Rq and brief phases above 
1OO7?0. The tidal radius of a star, equation (1), is a time- 
varying function that depends on the mean stellar density, as 

rt(t) oc M*(0"^/^^*(0 cx p*(0"^/^. The right panel of Figure 1 
illustrates the variation in stellar evolution profiles with initial 
mass and that these late phases of stellar evolution give rise to 
brief periods in which the star becomes orders of magnitude 
more vulnerable to tidal disruption than it is on the MS. 

2.3. Tidal disruption basics applied to evolved stars 

Tidal disruptions occur when stars are scattered into suf- 
ficiently low angular momentum orbits that they pass within 

1 /3 

the tidal radius at pericenter, r^ = (Mbh/M*) R^- In Figure 
2, we illustrate the respective tidal radii of stars in different 



evolutionary states. We denote the impact parameter of an 
encounter in terms of the ratio of the tidal radius, rt, to the 
pericenter distance, rp, as (3 = r^/r^, such that (3 I signi- 
fies a deep encounter. Tidal disruptions are only observable 
if they occur outside the black hole's Schwarzschild radius 
Ts = 2GMbh/c^. In terms of r^, the tidal radius may be rewrit- 
ten as 

For Mbh ^ 10^ Mq, solar- type stars will be swallowed whole, 
producing no observable flare. The precise value of this black 
hole mass cutoff is almost certainly modulated to some extent 
by the black hole's spin and innermost stable circular orbit 
(Kesden 2012; Haas et al. 2012). Although stellar remnants 
such as white dwarfs are expected to be numerous in galactic 
center environments (e.g. Alexander 2005), we do not con- 
sider them here since they will be swallowed whole by black 
holes Mbh ^ 10^ M0 (Luminet & Pichon 1989; Rosswog et al. 
2008a,b). 

A characteristic encounter timescale in tidal disruption 
events is the pericenter passage time. This timescale is equiv- 
alent to the stellar dynamical time, r^/vp ^ y/R^/GM = t^y^, 
for encounters with pericenter at the tidal radius, rp = rt, re- 
gardless of stellar structure or evolutionary state. Consid- 
ering encounters of varying impact parameter, the passage 
timescale becomes tp ~ /3~^^dyn- During the encounter, mate- 
rial is stripped from the stellar core and spread into two tidal 
tails. The material in one of the tails is unbound from the 
black hole and ejected on hyperbolic trajectories. The other 
is bound to the black hole and will return on a wide range 
of elliptical orbits. Following Rees (1988), we can derive 
the fallback time of the most bound material, ^fb, by assum- 
ing the star is initially on a parabolic orbit, £^orb = 0- Then, 

1 /3 

^min ^ GM^ /R^ (Mbh /M* ) , and the corresponding Ke- 
plerian period, ^fb = 27rGMbh(-2£^niin)~^^^, can be recast as 
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TABLE 1 

IAMq ZAMS initial models generated from mesa 



Fig. 2. — A cartoon, not drawn to scale, highlighting that tidal radius 
(shown as a dashed line) depends on the mean density of stellar material. 
As a star evolves off the main sequence and ascends the giant branch its tidal 
radius (and the corresponding timescales of a tidal encounter) increase by or- 
ders of magnitude. Shown here are representative MS (yellow), RG (orange), 
and AGB (red) stars being disrupted at respectively larger tidal radii. 



Accretion flares from tidaUy disrupted giant stars with ~ 
10- lOO/?0 are thus expected to rise and faU on timescales of 
years to hundreds of years. Since the total mass lost from the 
star is similar, the peak accretion rate M ~ AM/^fb for these 
events will be correspondingly lower than for stars on the MS. 

How often a star enters the zone of vulnerability delineated 
by the tidal radius depends on the stellar density and velocity 
dispersion. One can crudely estimate the rate of tidal disrup- 
tion as 



4/3 



100 km s" 



105 pc- 



lOORr. 



(6) 



by assuming that the black hole inhabits an isotropic sea of 
stars such that Ni^o ~ n^T^ta, and making use of the gravita- 
tional focus limit of the tidal disruption cross-section St (see 
Rees 1988, equation 2). In equation (6), the tidal disruption 
rate is linearly proportional to the tidal radius, Mso ^t- A 
star that spends most of its lifetime at IRq, 10% of its life- 
time at IORq, and 1% at IOOT^q would therefore have a sim- 
ilar likelihood of tidal disruption during each of these stages, 
even though their timescales span two orders of magnitude. 
However, as noted by Rees (1988), there are factors which 
complicate this picture. In particular, the isotropic descrip- 
tion is invalidated by the facts that the number density of 
stars falls of sharply with distance from the black hole and 
that stars within the black hole's sphere of influence orbit the 
black hole directly (in nearly Keplerian orbits). More detailed 
models find similar tidal disruption rates but different scalings 
than those derived under the simplifying isotropic hypothesis 
(Magorrian & Tremaine 1999; Wang & Merritt 2004). The 
implications of these more accurate models will be discussed 
in more detail in Section 4. 

3. HYDRODYNAMICS OF GIANT DISRUPTION 

The hydrodynamics of tidal disruption have been studied in 
detail for main sequence stars and planets, which are thought 
to be well described by poly tropes. Initial simulation ef- 
forts exploited semi-analytic formalisms that follow the dis- 
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Mcore 




^core 


^env 


/^core / A'-env 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


RGI 


1.4 


0.28 


12.5 


1.5 


1.7 


8.0 


0.048 


RGII 


1.32 


0.44 


110 


1.5 


1.8 


7.0 


0.053 


HB 


1.26 


0.55 


12.3 


1.5 


1.7 


8.0 


0.052 


AGB 


1.17 


0.57 


208 


1.5 


1.8 


7.0 


0.059 



Note. — Columns: (1) model name, (2) total mass [Mq], (3) core mass 
[Mq], (4) radius [Rq], (5) polytropic index of core, where F = 1 + 1/^, (6) 
polytropic index of envelope, (7) ratio of core to envelope mean molecular 
weights, (8) ratio of core radius to total radius. 

tortion of a sphere into a triaxial ellipsoid in the tidal field of 
the black hole (Carter & Luminet 1982; Luminet & Marck 
1985; Luminet & Carter 1986; Kosovichev & Novikov 1992; 
Ivanov & Novikov 2001). Later studies used both Lagrangian 
(Nolthenius & Katz 1982; Bicknell & Gingold 1983; Evans 
& Kochanek 1989; Laguna et al. 1993; Kobayashi et al. 2004; 
Brassart & Luminet 2008; Rosswog et al. 2009; Ramirez- 
Ruiz & Rosswog 2009; Lodato et al. 2009; Antonini et al. 
2011) and Eulerian formalisms (Khokhlov et al. 1993a,b; Di- 
ener et al. 1997; Guillochon et al. 2009, 2011; Guillochon & 
Ramirez-Ruiz 2012). In this section, we will emphasize some 
of the characteristics that make the disruption of giant stars 
unique. 

3.1. Initial models for hydrodynamic simulations 

We explore the typical characteristics of tidal disruption 
events involving a lO^M© black hole and stars in the repre- 
sentative stages outlined in Section 2 and shown in Figure 1 . 
To this end, we define 4 initial models based on the MESA 
evolution of a solar metallicity 1 .4M0 ZAMS mass star. We 
include two models for RG stars, one of which is ascending 
the RG branch (model RG I), and one at the tip of the RG 
branch (model RG II). A single model is used for each of the 
HB (model HB) and AGB (model AGB) phases of stellar evo- 
lution. 

For our hydrodynamic calculations we approximate the 
core and envelope structure using a nested polytropic structure 
with two fluids of different mean molecular weights (Chan- 
drasekhar 1967). This structure is approximate, but, as seen 
in Figure 3, the fits to the enclosed mass profiles within the en- 
velopes are quite good. We decided to enlarge the core radius 
relative to its true size for numerical convenience because it 
would be very restrictive to attempt to simultaneously resolve 
the physical and temporal scales of the true core (Rcoxt ^10^ 
cm and ^dyn ^ 5 s) while also following the evolution of the 
extended envelope (7?env ^ 10^^ cm and ^dyn ^10^ s). In our 
models we choose the core radius to be 7?core ^ ^env/20. This 
approximation imposes the restriction that we may consider 
only encounters for which the pericenter distance, rp, satis- 
fies rt(7?core,^core) < ^ n(^env,Menv). Thcsc are cncountcrs 
that are disruptive for the envelope while leaving the enlarged 
core tidally unaltered. In Table 1, we give the characteristics 
of our polytropic approximations to MESA models. 

3.2. Methods 

The calculations presented here have been performed with 
the tidal disruption code described in detail in Guillochon 
et al. (2009), Guiflochon et al. (2011) and Guillochon & 
Ramirez-Ruiz (2012). Our formalism is based on the 
FLASH4 code (Fryxell et al. 2000), an adaptive mesh Eule- 
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Fig. 3. — In dashed black, we show profiles of enclosed mass calculated in MESA for a IAMq ZAMS mass star at characteristic stages along its post-MS 
evolution. From left to right, these are: ascending the RG branch (RG I), the tip of the RG branch (RG II), HB, and AGB. The orange and blue lines show how 
we approximate the stellar structure using a nested, two-fluid polytrope. We match the profile of the envelope as accurately as possible while allowing the core to 
be enlarged Rcore ~ /20. Unless the passage is so close as to be disruptive for the core (/3 > 20), this approximation of the stellar structure is valid and lessens 
the computational burden of our calculations. 



rian hydrodynamics code. The compressibiHty of the gas is 
described with a gamma-law equation of state P oc . For the 
envelope we take j = j^^ = 5/3, while we choose a stiff equa- 
tion of state to describe the artificial core, 7 = 5. Because the 
core is only weakly compressible, it does not expand signifi- 
cantly as mass is removed from the star's envelope, ensuring 
that it makes no spurious impact on the dynamics of the tidal 
disruption. The flux between grid cells is solved using piece- 
wise paraboHc method (PPM) (Colella & Woodward 1984; 
Woodward & Colella 1984), which incorporates second-order 
directional operator splitting (Strang 1968). We base the grid 
mesh refinement on the value of the density, with successive 
cutoffs limiting the refinement by one level per decade in den- 
sity. 

Our calculations are performed in the frame of the star to 
avoid introducing artificial diffusivity (for a fixed resolution) 
by moving the star rapidly across the grid structure (Guillo- 
chon et al. 2009, 2011). The self-gravity of the star is com- 
puted using a multipole expansion about the center of mass 
of the star with /^ax = 10. The orbit is then evolved based on 
the center of mass of the star and the position of the (Newto- 
nian) point mass black hole. We refer the reader to Appendix 
A of Guillochon et al. (2011) for details of this algorithm. In 
our simulations the approximation of Newtonian gravity for 
the black hole is justified by the fact that the closest approach 
which any of the giant stars in our simulations make is about 
1 lOrs, well into the weak field regime. 

The simulations presented here are all resolved by at least 
7?^/Arnun > 200, where Armin is the dimension of the small- 
est cells. However, by implementing density cutoffs in the 
peak refinement for a given block we allow the linear reso- 
lution on the limb of the star to be 2 or 3 levels lower than 
the peak resolution of the core region. This limits our total 
discretization of the computational domain to (1 - 10) x 10^ 
cells, constraining the computational intensity of the problem 
while allowing us to survey some degree of impact parameter 
and stellar evolution parameter space. 

3.3. Pericenter passage and mass removal 

The passage through pericenter is mainly characterized by 
the impact parameter p = ri/r^. Since stars are generally scat- 
tered to the black hole from apocenter distances around the 
SMBH's sphere of influence (see Section 4), typical orbits 
leading to disruption follow nearly radial trajectories. Thus, 



we consider orbits which are initially parabolic. 

As the star passes through pericenter, quadrupole distor- 
tions of the stellar surface reach an amplitude of order unity. 
Material is sheared from the stellar envelope as the time vary- 
ing tidal field applies a gravitational torque to the distorted 
star, spinning stellar material to nearly its corotational angu- 
lar velocity. A portion of this spun-up material is unbound 
from the star and ejected into two tidal tails. The star's cen- 
ter of mass continues to follow a roughly parabolic trajectory, 
material in one of the tails is unbound and thus ejected onto 
hyperbolic orbits, and material in the other tail traces out ellip- 
tical orbits and will eventually fall back to the black hole. The 
tidal field of a lO^M© black hole is only ^3% asymmetric on 
the scale of a star at the tidal radius (the degree of asymme- 
try scales linearly with p). This symmetry is reflected in the 
formation of the two similarly shaped tidal tails as shown in 
Figure 4. 

While the structures and tidal deformations seen in Figure 
4 appear broadly consistent with those resulting from the dis- 
ruption of MS stars (Guillochon & Ramirez-Ruiz 2012), the 
quantitative nature of the distortions in the giant star case is 
marked by the gravitational influence of the dense core. Of 
course, MS stars have cores that are more dense than their 
envelopes, but the ratio of their central to average density is 
Pcore/p* ^ 100. whereas for giant stars this ratio is > 10^. In 
the initial models we construct, the density contrast is not as 
high due to numerical limitations, Pcore/ p* ~ 2 x 10^, but it is 
similarly marked by a discontinuous transition in the density 
profile. Because the dynamical time of the core is much less 
than that of the envelope, the core material is not perturbed, 
but its gravitational influence significantly modifies the rear- 
rangement of the surrounding envelope material as the en- 
counter progresses. The adiabatic response of centrally con- 
densed, nested poly tropes, like the ones constructed here, is 
to contract when they lose sufficient mass, thus shielding the 
star from further mass loss (Hjellming & Webbink 1987). The 
ability of the stellar core to retain envelope gas can be seen in 
the upper panels in Figure 4. These panels show a close view 
of the region near the star's core that leaves the passage tidally 
excited, but undisrupted. The core itself is completely unper- 
turbed. 

In Figure 5, we explore the causes of the mass loss in 
greater quantitative detail. The left panel of Figure 5 displays 
the star's self-bound mass as a function of time for a set of en- 
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Fig. 5. — Mass loss with varying impact parameter, /S, for the RG I model. The left panel shows the star's self-bound mass as a function of time, where t = 
corresponds to pericenter. The center panel shows the net mass loss, as measured after ~ SOt^y^, as a fraction of the envelope mass, Menv As /5 increases, the mass 
lost asymptotes to a value less than Menv- This can be contrasted to MS stars, which are fully disrupted at /5 ~ 1 .9. The reason for this asymptotic behavior can be 
seen in the adiabatic response of the stellar envelope in the expanding background of the core-black hole potential. In the right panel, we plot the enclosed mass 
within the (time-varying) Hill radius of the core, equation (7). The increase in enclosed mass after pericenter shows that the envelope is contracting, becoming 
more dense, and effectively shielding itself from further mass loss. 



counters with 0.75 < /3 < 2.25. Encounters of varying /3 lead 
to different degrees of mass loss and varying speed of mass re- 
moval (characterized by the slope of the lines near pericenter, 
^ = 0). The self-bound mass gradually decreases over many 
dynamical times, eventually converging to a net mass loss, 
-AM, plotted as a fraction of the envelope mass Menv in the 
center panel of Figure 5. The mass lost during the encounter 
asymptotes to less than the total envelope mass, with deeper 
encounters not being able to remove significantly more mass. 

The reason for the observed asymptotic behavior of -AM 
with (3 can be inferred from the right panel of Figure 5. The 
tidal force at pericenter increases with increasing (3, removing 
a greater fraction of the envelope mass as (3 increases. This 
increase in the rate of mass loss is drastically halted when 
the total envelope mass remaining approaches the mass of the 
core. Because the core is a significant fraction of the total stel- 
lar mass, the star's structure responds to a drastic loss of enve- 
lope mass by contracting. As more mass is removed, the sur- 
rounding envelope contracts more dramatically (Hjellming & 
Webbink 1987; Passy et al. 2011). The competition between 
the tidal mass stripping and the ensuing envelope contraction, 
both of which are seen to increase with f3, leads to the asymp- 
totic behavior with (3 of -AM. A subtlety lies in the fact that 



this contraction occurs not in an absolute sense, but in an ex- 
panding frame defined by the fluid trajectories in the com- 
bined potential of the black hole and the surviving core. This 
can be seen directly by measuring the enclosed mass within 
the Hill radius of the core, shown in the right panel of Fig- 
ure 5. The Hill radius of the core varies as function of time 
through the encounter as 

/?h,core(0= (^] ' r(0, (7) 

where r(t) is the separation between the core and the black 
hole. Menc(< ^h,core) decreases initially and reaches a mini- 
mum at pericenter (t = 0) where 7?h,core = /3~H^core/^*)^^^^*- 
After pericenter, the mass within this equipotential surface in- 
creases because envelope material contracts relative to the ex- 
panding frame of the disrupted star. 

Each of the post-MS models described in Section 3.1 dif- 
fers in the exact polytropic index of the envelope and the core 
to envelope mass ratio. However, the qualitative nature of 
mass removal from these models is very similar to that shown 
in detail for the RG I model in Figure 5. As can be seen in 
Figure 6, evolved stars retain an increasing fraction of enve- 
lope material at a given /3 as their core to envelope mass ratio 
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Fig. 6. — Mass loss as a function of /5 for the four post-MS models shown 
in Figure 3. Mass lost is shown as a fraction of the original envelope mass 
of a given star. Increasing core to envelope mass ratio dictates that that a 
smaller fraction of envelope mass is lost at a given impact parameter. How- 
ever, in all cases, much less than the entirety of the envelope mass is lost, 
even in encounters for which /S > 1. The solid lines are described by the 
fitting formulae presented in equation (8). 

grows. This is because larger core mass fractions in evolved 
stars increase the strength of the envelope's adiabatic response 
to contract with increasing mass loss (Hjellming & Webbink 
1987; Passy et al. 2011). In all cases, a significant fraction of 
envelope mass remains bound to the core for even the deepest 
encounters we have explored. This stands in stark contrast to 
MS stars (Guillochon & Ramirez-Ruiz 2012), which are fully 
disrupted beyond /3 ~ 1 .9. We tested both more and less dense 
artificial cores and confirmed that the asymptotic behavior of 
the mass loss is in no way affected by our enlarged core ap- 
proximation at the impact parameters we consider. 

The following fitting formulae provide a reasonable approx- 
imation to the degree of mass loss within the range of /3 pre- 
sented here 



-AM= < 



^ 0.52 + 1 . 1 1/3-^ - 1 .95/3-2 _^ 0.66/3-^ RG I, 

0.57 + 0.61/3-1 - 1.58/3-2 + 0.62/3-^ RG II, 

0.62 - 0.06/3-1 _ 0.30/3-2 _ o.04/3-^ HB, 

^ 0.55 - 0.15/3-1 - 0.32/3-2 + 0.01 fS'^ AGB . 



(8) 



Our findings may be compared to those of authors who 
study giant star stellar collisions (e.g. Bailey & Davies 1999; 
Dale et al. 2009). These studies similarly find that the dense 
core of the giant star retains a significant fraction of envelope 
mass. Dale et al. (2009) also note that the great majority of 
the star's hydrogen envelope must be removed to apprecia- 
bly alter the star's evolution along the giant branch. Because 
they retain sufficient portions of their envelope mass, the rem- 
nants of the tidal disruption events discussed here do not ap- 
pear as naked helium cores and instead are likely more similar 
in structure and appearance to their original, giant-star nature. 

3.4. Evolution of unbound material 

Material stripped from the surviving core forms two sym- 
metric tidal tails (see Figure 4). The evolution of the material 
in the tails is well described by two characteristic stages: an 
initial period of hydrodynamic evolution, followed by a pe- 
riod of homologous expansion. Hydrodynamic interactions 
remain important when the local sound speed of the stellar 
material exceeds the relative expansion velocity of adjacent 
fluid elements. As the stellar material expands, p decreases 



and the sound speed drops because = y^^yPjp ex p^^ 1^/2 

for 7 = 5/3. When these sound waves are unable to prop- 
agate as quickly as the distance between neighboring fluid 
parcels increases, the material's velocity distribution begins to 
freeze. Any additional evolution of the stellar fluid is then rel- 
atively well described by collisionless trajectories in the black 
hole's gravitational field. 

To determine the trajectories of the material in the tails, we 
calculate the evolution of the specific binding energies of all 
of the fluid elements relative to the black hole. This gives 
the spread of mass per unit energy within the disrupted star, 
dM/dE, where E is the specific orbital binding energy rela- 
tive to the black hole. In Figure 7, we plot the time evolu- 
tion of dM/dE with the self-bound material of the surviving 
stellar core excised. The majority of the material achieves ho- 
mologous expansion only after a few dynamical timescales, 
with the structure of the material near the tips of the tidal tails 
being frozen in earliest. The distribution of material close 
to the the surviving core continues to evolve for some time 
(center panel of Figure 7); it will not fully converge until the 
surviving remnant is truly isolated from the black hole. This 
effect can also be seen in the left panel of Figure 5, where the 
self-bound mass of the surviving star decreases gradually as 
a function of time. The gentle evolution of the stream ma- 
terial at the force balance point between the remnant and the 
black hole occurs over a long timescale when compared to the 
pericenter passage time because the surviving star continues 
to be slowly depleted of mass that crosses the lowered poten- 
tial barrier between the star and the black hole. Hence, the 
time evolving gravitational influence of the surviving star has 
to be taken into account to accurately compute the resulting 
dM/dE distribution. 

The material in the tails follows ballistic trajectories after 
sufficient time has elapsed. As a result, the spread in energy 
can be mapped to a return time to pericenter for the tail that is 
bound to the black hole. 



M = 



dM dMdE (27rGMbh)2/^ 



dt dE dt 



dM 
~dE 



-5/3 



(9) 



For a flat dM/dE, one recovers the canonical r^^^ power law 
predicted for the fallback mass return rate after a tidal disrup- 
tion event (Rees 1988; Lodato et al. 2009). The right panel 
in Figure 7 shows the mapping to M from dM/dE plotted in 
the left two panels. The cavity in dM/dE that results from 
not including self-bound material leads to a similar cavity in 
the predicted fallback rate at very late times. As the rem- 
nant flies away from the black hole, becoming increasingly 
isolated from the black hole's potential, we expect this distri- 
bution to fill in (but not necessarily completely flatten). 

3.5. Debris fallback and AGN flaring 

Debris ejected into the tidal tail bound to the black hole 
returns to pericenter at a rate given by equation (9). The re- 
turning gas does not immediately produce a flare of activity 
from the black hole. First material must enter quasi-circular 
orbits and form an accretion disk (Ramirez-Ruiz & Rosswog 
2009). Once formed, the disk will evolve under the influ- 
ence of viscosity (Cannizzo et al. 1990; Montesinos Armijo 
& de Freitas Pacheco 2011). However, the viscosity would 
have to be extremely low for the bulk of the mass to be stored 
for longer than (equation 5) in a reservoir at r ?^ rt. Taking 
typical values, the ratio of the viscous accretion timescale to 
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-2x1015-1x1015 1x10^5 2x10^5 
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Fig. 7. — Temporal evolution of the specific binding energy of the tidal tail material, calculated here for the AGB model at an impact parameter of /3 = 1. The 
lines plotted range from t = St^yn (blue) after pericenter tot = 35fdyn (brown) after pericenter in intervals of St^yn and include only the material which is unbound 
from the surviving core. In the left and center panels, specific binding energy relative to the black hole is shown on on the x-axis. The star approaches the black 
hole on a parabolic orbit corresponding to £^ = 0. Material with £^ > is ejected on hyperbolic orbits, while that with £^ < eventually falls back to the black hole 
on a range of elliptical orbits that may be mapped to the projected fallback rate shown in the right panel using equation (9). The energy distribution of material 
near the tips of the tidal tails is frozen into homologous expansion at early times. Material closer to £^ ~ evolves to fill in the cavity in dM/dE for many r^yn 
as the surviving remnant becomes increasingly isolated from the black hole's potential. The dM/dE and M distributions shown in Figures 7 - 9 are plotted with 
500 bins in energy or time. 




Fig. 8. — The fallback rate of bound material with varying (3 for 
the RG I model, normalized to the Eddington accretion rate, Medd = 
O.O2e-\(Mbh/lO6M0)M0yr-i. This may be compared to Figure 5, which 
details the removal of mass through pericenter for this model. Differing im- 
pact parameter changes Mpeak and fpeak to some extent. In deeper encounters, 
the material most bound to the black hole falls back earlier and exhibits a 
more gradual rise to peak. After peak, all of the /3's shown here appear to fall 
off more steeply than t~^/^ for many years. The other three post-MS models 
exhibit very similar qualitative behavior. 
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Fig. 9. — A comparison of the fallback accretion between the four evolved 
star models and the disruption of a sun-like MS star, at /3 = 1.5. Accretion 
flares from the disruption of giant stars are long-lived and peak at relatively 
lower M. While the 100-1000 year typical timescales of stars disrupted at the 
tip of the giant branches will be difficult to distinguish from other long-term 
low-level AGN activity. The rise to peak resulting from the disruption of stars 
either ascending the red giant branch (RG I) or on the horizontal branch stars 
(HB) should be observable in current time domain surveys. We extrapolate 
the fallback beyond the position of the cavity in M (dashed lines) using the 
slope of M captured by our last simulation snapshot. 



the fallback timescale is 

% ^ 1 10-2 J [k^MqJ [iMeJ ' 

where we have made the standard assumption of a thick 
(H/R) ~ 1 disk (Shakura & Sunyaev 1973). The accretion 
rate onto the black hole, and to a certain extent the bolometric 
luminosity of the resulting flare (i.e. L ex Mc^), is thus gen- 
erally expected^ to be limited by the gas supply at pericenter 
M, equation (9), not by the rate at which the orbiting debris 
drains onto the black hole (Rees 1988; Ulmer 1999). 

In the simulations of giant star tidal disruption presented 
here, we find sizable variations in the black hole feeding rate, 
M, with impact parameter (3. In Figure 8, we show M as a 
function of /3 for the RG I model. The time of peak, ^peak, is 
observed to vary only slightly with (3. Additionally, since AM 

^ When the black hole is fed above the Eddington limit (M > Medd). the 
relationship between M and disk luminosity is, however, less certain (e.g. 
Strubbe & Quataert 2009, 2011). 



asymptotes at high /3, the peak accretion rates, Mpeak, are sim- 
ilar for most deep encounters. The slope of the initial fallback 
after peak is, however, significantly steeper than r^^^, par- 
ticularly for the more grazing encounters. This variation in 
post-peak slope differentiates M curves that overlap at peak. 

The time of initial fallback, ^fb, varies considerably with /3. 
Because the time of peak, ^peak, is relatively constant, the slope 
of the early-time rise varies correspondingly, as can be seen 
in Figure 8. Observations of this portion of the M curve could 
help break the relative degeneracy seen in the post-peak light 
curves. However, at these early times viscous and angular 
momentum redistribution processes might play some role in 
shaping the accretion luminosity. Further, because the early 
portion of M represents only a small fraction of the total mass 
lost by the star, AM, it is sensitive in our simulations to the 
limited mass resolution of discretized material in the tidal 
tails. In Figure 7, we show that the early-time M freezes into 
homologous expansion quickly, and is seen to accurately pre- 
serve its shape as the disruption ensues. We have studied the 
variations in the early time M with increasing resolution and 
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logio(Mbh) [Mq] 

Fig. 10. — Encounters of /5 = 1 rescaled to the range of known super- 
massive black hole masses. With increasing black hole mass, flares become 
longer-lived but also peak lower relative to the Eddington rate for that black 
hole mass. The forbidden region shows where rt < n and stars enter the black 
hole's horizon whole without being tidally disrupted. 

find no systematic variations. We do note, however, that the 
short-timescale bumps and variations in M at ^ ~ ^fb Hkely re- 
sult from grid (de-)refinement events, which are unavoidable 
in keeping our computations feasible. 

While variations in M with (3 are present, major changes are 
observed in the feeding rates as stars evolve and their struc- 
tures are dramatically altered. In Figure 9, we compare the 
black hole feeding rates from the disruptions of evolved stars 
with a typical MS disruption. The disruptions of giant stars 
lead to long-lived flaring events that peak at, or near, the Ed- 
dington rate for a lO^M© black hole. The peak time of these 
events is, as expected, correlated with the radial extent of the 
evolving star. Our simulations show that while the timescale 
of peak accretion rate in RG I and HB star flares is several 
years, the rise toward peak M occurs over similar timescales 
to those seen in the disruption of sun-like stars, thus offering 
hope for detectability in surveys tailored to detect MS flares. 
On the other hand, tidal disruptions of stars at the tips of the 
giant branches (the RG II and AGB models) lead to flares that 
peak at hundreds of years timescales, which make them diffi- 
cult to discern from other non-transient AGN feeding mech- 
anisms. In Figure 10, we show how the peak M and fallback 
timescales are expected to scale with increasing black hole 
mass. 

4. LOSS CONE THEORY AND RATES OF GIANT DISRUPTION 

Following the observation by Hills (1975) that stars passing 
within rt at pericenter are tidally disrupted, a series of studies 
focused on the rate at which stars are fed into disruptive orbits. 
This rate depends on the distribution of stars around the black 
hole and the rate at which small scatterings in orbital angular 
momentum allow stars to diffuse into nearly radial loss cone 
orbits that lead to disruption (Frank & Rees 1976; Lightman & 
Shapiro 1977). Cohn & Kulsrud (1978) present a formalism 
for treating this diffusion process by numerical integration of 
a Fokker-Plank equation. More recent work extends the Cohn 
& Kulsrud (1978) formalism to compute the flux of stars into 
the loss cones of observed galaxies (Magorrian et al. 1998; 
Magorrian & Tremaine 1999; Wang & Merritt 2004), with a 
key result being that stars are typically fed into the loss cone 
from weakly bound orbits with semi-major axes similar to the 



black hole's sphere of influence, equation (11). Magorrian 
& Tremaine (1999) also considered the effect of an axisym- 
metric (rather than spherical) stellar distribution, while Wang 
& Merritt (2004) considered the effects of black hole mass 
on tidal disruption rates. Additional recent work has focused 
on the effects of black hole spin (Kesden 2012), binary black 
holes (Ivanov et al. 2005; Chen et al. 2009), and numerical 
determination of the tidal disruption rate using n-body calcu- 
lations (Freitag & Benz 2002; Brockamp et al. 2011). 

Syer & Ulmer (1999) use a simplified treatment of both the 
loss cone formalism and stellar evolution to roughly account 
for the effects of evolved stars on the integrated disruption 
rate. We extend the work of previous authors by emphasizing 
the relative rates of disruption of stars in different evolution- 
ary phases using our MESA stellar evolution models, but we 
otherwise follow the Fokker-Plank formalism closely (Cohn 
& Kulsrud 1978; Magorrian & Tremaine 1999; Wang & Mer- 
ritt 2004). In this section, we begin by describing a simplified 
galactic nuclear cluster model. We outline the loss cone for- 
malism and compute how the tidal disruption rate changes as 
a star's tidal radius changes due to stellar evolution. We com- 
pare our findings based on the simple nuclear cluster model to 
a sample of observed nuclear cluster profiles to infer the ef- 
fects of structural differences between galaxies. We are then 
able to use our calculation of how the tidal disruption rate 
scales with tidal radius to determine the stellar ingestion diet 
of SMBHs. 

4. 1 . Simplified nuclear cluster model 

To develop intuition for the relative rates of disruption of 
MS and evolved stars, we first explore a simplified nuclear 
cluster model. This model consists of a homogeneous stel- 
lar population and is described by a Keplerian potential. The 
black hole is the dominant influence on the stellar kinematics 
within a nuclear cluster of radius 



lO^Mrr 



1/2.12 



pc, 



(11) 



where is the external velocity dispersion of the greater 
galactic bulge. In the numerical expression above we use 
the Mbh - cr relation (e.g. Ferrarese & Merritt 2000; Geb- 
hardt et al. 2000; Tremaine et al. 2002; Gultekin et al. 2009), 
with fitting values from Gultekin et al. (2009), = 2.43 x 
1O5(M0/Mbh)'^-^'^ cm s-\ 

At radii r < rh, the cluster potential is approximately Kep- 
lerian, 

^(r)= , (12) 

r 

where, following Magorrian & Tremaine (1999), we choose 
a positive sign convention for the potential and the energies 
of bound orbits. The velocity dispersion at a given radius is 
approximately cr^(r) ~ ^(r). We define stellar orbits based on 
their binding energy to the black hole 



£(a) = 



GM\ 



bh 



la 



(a), 



(13) 



where a is the orbital semi-major axis. The Keplerian orbital 
period of stars with a given e is P{e) = 27rGM\yh(2e)~^^^. 

We take a singular isothermal sphere stellar number density 
profile. 



y*{r) = u^r 



(14) 
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where = ^*(^h), as Wang & Merritt (2004) do. This is com- 
parable to the steep density cusp expected to form in a relaxed 
stellar distribution about a black hole (Bahcall & Wolf 1976, 
1977), while also being consistent with the roughly flat veloc- 
ity dispersion profiles observed in typical galaxies outside r^. 
For simplicity, a single stellar mass is commonly assumed, 
and, as a result the mass density is related to the number den- 
sity by p*(r) = M^v{r). We choose so that the mass of stars 
within the sphere of influence is twice the black hole mass, 
which gives = (27rr^)~^(Mbh/M^). We further assume that 
stellar orbits are isotropic in angular momentum space and 
therefore may be described by a distribution function solely 
in energy. 



2x-3/2 ^(3) 

r(i.5) 



1/2 



(15) 



where the numerical factor r(3)/r(1.5) ^ 2.26 (equation 9 in 
Magorrian & Tremaine 1999). Having described the proper- 
ties and phase space distribution of stars in our cluster model, 
we can now determine how often stars are disrupted. 

4.1.1. The full loss cone disruption rate 

The angular momentum of a circular orbit is the maximum 
angular momentum allowed for a given e. It is 



The angular momentum of a loss-cone orbit is 



y2(£) = 2r; 



.2 

min 



GMi 



bh 



' 2GMbhr^ 



(16) 



(17) 



where rmin is the maximum of or r^. 

The number of stars in a full loss cone per energy is then 
Nicie) = 47r^ f(e)P(e)Jl(e), and the total number is 



Mc = f 

Jo 



(18) 



where Smax is the energy corresponding to rmin. The flux of 
stars into the loss cone, if we assume that it is continuously 
replenished, is simply 



FMi(e) = Ni,(e)/P(e). 



(19) 



Similarly, the net rate at which stars are disrupted from a full 
loss cone is obtained by integrating over energy. 



(20) 



However, the loss cone is not necessarily full for all e. We 
turn or attention to where Nf^w is an appropriate estimate of 
the tidal disruption rate. 

4.1.2. Repopulating loss cone orbits 

The derivation of Nf^w assumes that at all energies, £, stars 
destroyed by the black hole at pericenter are replaced on a 
timescale shorter than the orbital period, P{e), thus ensuring 
that the phase space density of low / orbits remains unde- 
pleted. A typical timescale for stars to random walk in an- 
gular momentum by J\c via the two-body relaxation process 
is 



hie) ' 



Me) 



(21) 



where two-body relaxation time is 

5.68 X 10^9 Mo lO^Mo pc"^ 



tr(r) = 



In A p*(r) 



^(r) 
lO^cm/s 



(22) 



(equation 7-107 in Binney & Tremaine 2008). The Coulomb 
logarithm is approximated as A ^ 0.4A/'*(r < rh), where 
A/'*(r < Th) is the number of stars within rh, or equivalently, 
the number of stars with e > e^- 

Thus, in an orbital period, the typical change in angular mo- 
mentum a star receives via this random walk process is 



Aj\e) = Jl(e)P(e)/tj(e), 



(23) 



or, written as a ratio to the loss cone angular momentum 
(Cohn & Kulsrud 1978), 



(24) 



We can imagine two limiting cases for the value of q. The 
full loss cone rate is appropriate when AJ J\c. This is the 
so called pinhole limit where per-orbit scatterings are large 
compared to the target, here given by the size of the loss cone 
(Lightman & Shapiro 1977). Loss cone orbits are easily re- 
populated in an orbital time in this limit, ensuring that Fmi 
is satisfied. In the opposite limit, AJ <C J\c, the loss cone 
refills much more gradually than it is depleted. Once in the 
loss cone, stars have little chance of being scattered out in an 
orbital period. They are disrupted as they pass through peri- 
center and the loss cone remains mostly empty. This is known 
as the diffusion limit, where stars must diffuse toward the loss 
cone over many orbital periods (Lightman & Shapiro 1977). 

The flux into the loss cone in the diffusion limit is less than 
the full loss cone rate, since the loss cone will be depleted. We 
define the parameter Rq to specify the low angular momentum 
limit above which orbits are typically populated. 



Ro(e) 
Ride) 



^ / exp(-^) 
\ exp(-0.186^-0.824y^) 



if^(^:)> 1, 
ifq(s)<l. 



(25) 



where R = J^/J^ (Magorrian & Tremaine 1999). This allows 
us to calculate the loss cone flux in the transition between 
the two limiting cases, where loss cone orbits are partially 
repopulated in an orbital period. The flux of stars into the loss 
cone is given by 



Fide) = 47r^AjHe) 



m 

HR-,') 



(26) 



(Magorrian & Tremaine 1999). In the pinhole limit, Ro(s) 
0, and F\c(£) ^fuii(^)- For more tightly bound orbits corre- 
sponding to the diffusion limit we have Roie) Rid^) and 
^ic(^) ^fuii(^)- Finally, the integrated flux of stars into the 
loss cone is obtained by integrating F\c(£) over e, 



N = 



F[c(e)de. 



(27) 



Figure 1 1 shows the loss cone flux for the simplified nuclear 
cluster model described here. The peak of the flux comes 
from e ^ Eh, which corresponds to orbits with semi-major 
axes comparable to the black hole's sphere of influence, rh. 
At energies lower than the peak, the pinhole limit is applica- 
ble, and F\c ~ Fm\. At high energies, F <C i^fuii because the 
loss cone is only partially refilled each orbital period. 
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Fig. 11. — Loss cone flux Fic of solar-type stars for the simplified nuclear 
cluster model of Section 4.1. Orbital binding energy is inversely proportional 
to semi-major axis £ oc I /a, and £h is the energy corresponding to the radius 
of influence of the black hole, in this case rh ~ 1 parsec. At low e, AJ ^ Jic, 
and the loss cone is refilled on an orbital timescale. In this limit, F\c follows 
the full loss cone rate. At higher e, AJ <C J\c, and orbits with / < J\c are 
incompletely repopulated in an orbital period. In this limit, F[c is set by the 
diffusion rate in orbital angular momentum. At £ ^ £h a star must complete 
A/orb ~ /ic/A/ orbits to diffuse one loss cone angular momentum, J\c. At each 
pericenter passage stars are vulnerable to collisions, leading to a collisional 
truncation of F[c where A/orb ^ 1 • 



4.1.3. Modification ofF\c by direct stellar collisions 

In the crowded nuclear cluster environment, the large ge- 
ometric cross-section of giant stars leaves them susceptible 
to collisions with other stars. The effect of these collisions on 
F\q is most pronounced in the diffusion limit, where stars must 
complete A^orb ^ - orbits to be able to diffuse 

a total A/ = /ic in angular momentum space. To estimate the 
importance of collisions, we compute the probability of sur- 
vival, Psurv, during one angular momentum diffusion time, t]. 
The collision cross section is 



1 + 



2G(Mi+M2)' 



(28) 



where Rj^m = /^i +^2 and the two stars are represented by sub- 
scripts 1 and 2. We approximate low angular momentum or- 
bits (which may diffuse to loss cone orbits within ^j) as radial 
paths, s, through the cluster from apocenter to pericenter and 
back again. The typical number of collisions suffered during 
tj is then 



A^coll(^j) = 



q(e) 



n2a{E) 

L ' 



j:(Ri-hR2,sMs)ds, 



(29) 



and the probability of survival is /^surv(^j) = 
Max[0, 1 -7^A/coii(^j)] , where r] represents the efficiency 
of collisions in destroying stars. For our purposes, 7^ = 1 is 
a reasonable choice because even if stars are not completely 
destroyed by a collision, they suffer a large angle scattering 
sufficient to deviate their orbital parameters significantly 
from their previous random walk. 

The importance of collisions can be seen in Figure 11 
(which considers solar- type stars). At high orbital binding 
energies, e, F\c is truncated by the strong likelihood of stars in 
low / orbits suffering a collision in t]. 

4.1.4. Scaling of the loss cone flux with tidal radius 

To determine the relative contribution of different stellar 
evolutionary stages to the tidal disruption rate, we consider 
the per star (specific) tidal disruption rate ti. The specific rate, 
ri, varies over the stellar lifetime with changing tidal radius. 



and, because we consider a cluster of homogeneous stellar 
constituents, n = N/N. In Section 2.3, we show that for an 
isotropic distribution of stars, the flux into the loss cone scales 
linearly with the tidal radius. This is only true because the loss 
cone is full, thus the rate depends on the size of the loss cone, 
and n (X J^^ (X r^. If the loss cone is empty, changes in its size 
do not have an effect on the rate, which is set instead by the 
diffusion process. In the terminology of Lightman & Shapiro 
(1977), these are the pinhole and diffusion limits, respectively. 
Taking the relevant limits of equation (26), we have 



ln(/2)aln(rt) 



for q(s) 1 (pinhole), 
for q(e) <C 1 (diffusion). 



(30) 



The exact scaling the tidal disruption rate, ri, must be some 
combination of the contribution from these two limiting cases 
because it is an integral quantity in e that spans both limits. 

The effect of the superposition of the full and empty loss 
cone limits within a particular nuclear cluster is shown in Fig- 
ure 12. The left panel shows the change in the loss cone flux, 
Fic, as Tt increases by a factor of 10 and 100. In the portion 
of the cluster where F\c ~ i^fuii, ^ic scales linearly with r^. In 
the empty loss cone limit F\c is approximately constant with 
increasing r^, as expected. Because Jic is increasing while 
AJ remains the same, the energy corresponding to q(e) = 1 
(roughly the peak in F\c) moves to lower e. Therefore, as the 
loss cone increases in size, stars are typically fed from less 
bound orbits that originate further from the black hole. 

Specific tidal disruption rates, ri, can be calculated via the 
integrated loss cone flux, N, equation (27). This integral is 
shown for a range of rt and Mbh in the right panel of Figure 
12, normalized to the specific rate of solar- type stars for a 
given Mbh. For a wide range of tidal radii, the integrated rate 
follows ri (X r^, with a = 1/4. At large black hole masses 
is larger than rt,©, and, as a result, solar- type stars are swal- 
lowed whole rather than producing a disruption, as described 
by equation (4). This is coded in the right panel of Figure 12 
as grey shaded regions for stars consumed by the black hole 
and blue shaded regions for flaring events. For stars larger 
than ~ 1OO7?0 collisions become a dominant effect and n sub- 
stantially decreases. 

4.2. Impact of the structural diversity of observed galactic 

centers 

In Section 4.1, we compute the relative likelihood of tidal 
disruption for stars of varying tidal radii under the assumption 
that the stellar density distribution is accurately described by 
a singular isothermal sphere. However, the central regions of 
observed galaxies exhibit somewhat more complex structure. 
Out of this complexity arises some degree of galaxy to galaxy 
variation, which we illustrate here. 

4.2.1. Galaxy sample 

Observations resolving the central SMBH's sphere of influ- 
ence in nearby, early type galaxies have been performed by 
the Nuker team using the Hubble Space Telescope (see Lauer 
et al. (1995) for early work and Lauer et al. (2007) for a more 
recent review). The radial surface brightness profile in these 
galaxies is found to be bimodal (Faber et al. 1997). Many 
of the most massive galaxies exhibit a core or flattening of 
the surface brightness profile at small radii. The presence of 
such cores is thought to be lingering evidence of gravitational 
heating from the inspiral of a pair of black holes following a 
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Fig. 12. — Scaling of the tidal disruption rate with increasing tidal radius. The left panel shows that with increasing n in the A/ ^ J\c limit, F\c scales linearly 
with rt, while in the A J <C J\c limit, F[c scales only weakly with n. We assume that increasing stellar radius is the mechanism causing rt to increase (because 
in stellar evolution the changes in mass are small). As a result, for larger rt collisions become increasingly important. The right panel shows the integrated rate 
per star, n. For a wide range of tidal radii, the scaling is well described by a power law with slope, a = 1/4. At very large rt collisions become dominant and n 
decreases, while for small rt stars are typically swallowed whole by the black hole because > rt, and n no longer depends on rt. The shaded regions describe 
the relative fractions of swallowed (grey) versus flaring (blue) events for a range of black hole masses between lO^M© and lO^^Mo. 



major galactic merger (Faber et al. 1997). The surface bright- 
ness profile of less massive galaxies typically continues to rise 
at smaller radii up to the resolution limit of the observations 
and is well described by a single power law. These galax- 
ies are termed power-law and exhibit central stellar densities 
tens to thousands of times higher than those inferred for core 
galaxies (Lauer et al. 2007). Nuker galaxy surface brightness 
profiles are fit to a parameterized, smooth function (the Nuker 
law) that describes the asymptotic slopes of the surface den- 
sity profile at large and small radii relative to a break radius 
Tb (Byun et al. 1996). The rates of tidal disruption events 
in these galaxies are computed by Magorrian et al. (1998) 
and Magorrian & Tremaine (1999) and revised by Wang & 
Merritt (2004) to accommodate changes in inferred black hole 
masses. 

4.2.2. Scaling ofn with n 

To investigate the influence of the bimodal nature of galac- 
tic center profiles on the scaling of fi with rt, we follow the for- 
malism of Magorrian & Tremaine (1999) and Wang & Mer- 
ritt (2004). We use the galaxy sample from Wang & Merritt 
(2004), who used data from Faber et al. (1997) and fit black 
hole masses based on the Mbh - cr relation of Merritt & Fer- 
rarese (2001). We show the deprojected mass density profiles 
of the galaxy sample used in our analysis in the left panel of 
Figure 13. At radii less that the break radius, rb, the mass 
density profiles of core and power-law galaxies separate dra- 
matically. Black holes with the highest inferred masses are 
typically found in core galaxies, rather than power-law galax- 
ies. The center panel of Figure 13 shows the scalings of the 
specific tidal disruption rate, fi, which is plotted scaled to the 
specific disruption rate of solar type stars Uq . Lines are col- 
ored here by the galaxies classification by Faber et al. (1997) 
as core or power-law, and we ignore collisions for simplicity. 
Where the tidal radius is less than the Schwarzschild radius, 
we still compute the ^(rt), but color the lines in gray. 

We find substantial variation between galaxies: a 1007?© 
star is between 2 and 10 times more likely to be disrupted than 
a solar- type star depending on the galaxy in which it resides. 
The dashed lines in the center panel of Figure 13 show power 
law scalings fi oc with a = 0,0.25,0.5,0.75, and 1. Many 



galaxies fall in the range of a ~ 0.2-0.5, broadly consis- 
tent with the simplified model presented in Section 4.1, which 
found a = 1/4. Another effect is a general trend toward shal- 
lower slopes dln(fi) / dln(rt) as rt increases, as seen in the right 
panel of Figure 13. The change in slope is most significant in 
core galaxies with massive black holes. The peak in loss cone 
flux moves out in the cluster as rt increases, feeding stars from 
regions of lower stellar number density. In core galaxies, the 
peak in Fic may cross rb with significant effect on the scaling 
of fi arising from the change in slope of the stellar number 
density at the break radius. In general, we find that the tidal 
disruption rate appears to vary more from galaxy to galaxy 
than between core and power-law classes. In what follows, 

we will use the relatively weak scaling fi ex rt as a value that 
is representative of the majority of galaxies, irrespective of 
their structure. 

4.3. The stellar diet ofSMBHs 

We have calculated how the specific tidal disruption rate fi 
changes as a function of rt. The tidal radius of a star is a func- 
tion of the initial stellar mass and age (see Figure 1), so one 
must integrate over a stellar lifetime to find the relative disrup- 
tion probability for stars during different evolutionary stages. 
Because ^ is a cluster-integrated rate, we are implicitly assum- 
ing that stars of different evolutionary stages are distributed 
isotropically in energy space. The expectation value for the 
number of disruption events a particular star will experience 
during an evolutionary stage lasting from ti to t2 is given by 
the integral of ^(rt) across that time period, 

nu = / fi{n)dt, (31) 

where rt is a function of M^{t) and R^{t). Similarly, the total 
lifetime integral is 

^tot= / fi{n)dt, (32) 

?ZAMS 

where ^zams is the age of the star at the zero-age main se- 
quence, and ^max is the age at the end of the star's lifetime. 
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Fig. 13. — Scaling of the tidal disruption rate in a sample of 41 early type galaxies whose central surface density profiles were observed with //^Srand fit with 
parameterized Nuker laws (for the details of the sample see Wang & Merritt 2004). In the left panel, we show the mass density of stars as a function of the 
density at the break radius, rb, a parameter of the two-power Nuker fit. Inside rb, galaxies show bimodal structures, with those hosting the most massive black 
holes typically having low density cores, whereas those with less massive black holes typically have dense central regions and single power law profiles. In the 
center panel, we plot the scaling of the specific tidal disruption rate n with n for each of these galaxies, normalized to the rate of solar type stars, uq . Where 
the tidal radius is less than the Schwarzschild radius we still compute the rate for illustrative purposes, but color the lines in gray rather than by galaxy structure. 
Finally, the right panel shows the slope of the scaling curves plotted in the center panel as a function of n. There is more substantial galaxy to galaxy variation 
in the slopes of these scaling curves than there is between structural class of galaxies. However, in most cases, <iln(^)/<iln(rt) <C 1 and in many cases is broadly 
consistent with the a = 1/4 scaling derived in the simplified model of section 4.1. 



Therefore, the fractional likelihood of disruption during a 
given evolutionary stage as a function of initial stellar mass 
is given by 

/stage(MzAMs) = • (33) 

^tot 

We can consider the probability of being disrupted at differ- 
ent evolutionary stages to be divided in this manner as long 
as most stars of a given Mzams do not experience a disruption 
during their lifetime (^tot < !)• The fact that we observe gi- 
ant stars in our own galactic center is a testament to the fact 
that many stars must survive without being disrupted during 
their entire lifetime (Genzel et al. 2010). This requirement is 
most likely to be satisfied for stars that are only weakly bound 
to the black hole (of which there are many and for which the 
black hole is an extremely small target). In this case, which 
corresponds to the peak in the feeding rate of giant stars, a 
typical star can easily survive without being disrupted over 
its entire lifetime, and we have N /N ^ ^max, indicating that 
^tot <^ 1- However, our assumption is less justified when con- 
sidering stars well inside the sphere of influence. Within this 
smaller reservoir of stars there is a relatively high probability 
of disruption integrated over the star's lifetime. In this case, 
there is a non-negligible chance that the star would have al- 
ready suffered a tidal interaction with the black hole before 
reaching a given evolutionary stage - either by evolving onto 
the loss cone (Syer & Ulmer 1999) or by the typical diffusion 
process. 

The fractional likelihood of tidal disruption during different 
evolutionary stages, equation 33, is shown in Figure 14 for a 
range of initial stellar masses, Mzams- We show how this flar- 
ing diet depends on the SMBH mass and the scaling of fi with 
rt. In calculating /stage, we assume that the stellar population 
is old enough that stars are able to reach each of these rep- 
resentative evolutionary stages. The top panels in Figure 14 
show the relative fractions of disrupted stars at different evo- 
lutionary stages given a ^ ex /J^ scaling, which we found to 
be representative of most galactic centers. The bottom pan- 

1 12 

els show the result given the optimistic scaling of ^ ex r/ . 



For black holes of relatively low mass (the sort that typically 
inhabit power-law galaxies), disrupted stars are typically MS 
stars with giant stars contributing 10-40% of all disruption 
events. As black hole mass increases, an increasing fraction 
of MS stars are consumed whole rather than disrupted. The 
flaring diet of the most massive black holes consists nearly 
exclusively of post-MS stars. The disruption fractions corre- 
sponding to a population of stars can be constructed using the 
estimates of /stage given in Figure 14 by weighting them with 
the stellar population mass spectrum and age distribution. 

5. DISCUSSION 

We have shown that the tidal disruption of stars while they 
are on the giant branch gives rise to longer-lasting accretion 
flares than when disruption occurs during the MS (Figure 9). 
Additionally, the stellar ingestion flaring diet of the most mas- 
sive SMBHs consists exclusively of evolved stars (Figure 14). 
The question remains whether these long duration accretion 
events occur frequently enough for flaring activity to be a dis- 
cernible characteristic of the most massive SMBHs. In this 
section, we use Monte Carlo realizations of flaring events to 
illustrate the defining properties and expected rates of flaring 
events as a function of black hole mass, their contribution to 
the accretion luminosity of local AGN, and the prospects for 
detecting AGN flares powered by the disruption of evolved 
stars. 

5.1. A limiting SMBH mass scale for tidal disruption flares 

Giant stars are the only stars that can light up the most mas- 
sive SMBHs by tidal disruption, since MS stars are swallowed 
whole rather than being tidally disrupted and subsequently 
producing a luminous accretion flare. Giant disruptions are, 
however, relatively infrequent due to the combination of the 
short lifetimes of evolved stars and the weak scaling of the 

tidal disruption rate with tidal radius fi cx /J^ . Therefore, de- 
spite their large radii, evolved stars are unable to compensate 
for the drop in flaring rate at black hole masses where MS 
stars are swallowed whole. As a result, we find a sharp cutoff 
in the tidal disruption flaring rate around Mbh ^ lO^M©. This 
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Fig. 14. — The stellar ingestion flaring diet of SMBHs is shown for a range of SMBH masses, and two scalings of n with rt, n (x . The fractional likelihood 
of tidal disruption during different stellar evolutionary stages, /stage, is plotted as a function of initial stellar mass, Mzams- Stars more massive than about 1 
solar mass evolve to the post-MS in less than a Hubble time and these evolved stars contribute a significant fraction of disruption events. As Mbh increases, an 
increasing fraction of MS stars are swallowed whole rather than being tidally disrupted. Flaring events on the most massive black holes consist exclusively of 
post-MS stars. 
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Fig. 15. — In the left panel, we show the fractional composition of stars scattered into the loss cone, f^, for a range of black hole mass between lO^M© and 
1O^^M0. We show a sharp cutoff in the flaring rate for black holes more massive thanlO^M©, beyond which MS stars are swallowed whole rather than tidally 
disrupted. Post-MS stars contribute some flaring events above this cutoff, but the weak scaling of n with tidal radius (Figure 12) dictates that these events are few 
compared to the flaring rate of lower mass black holes. The right panel shows the demographics of the flaring fraction only. These plots are based on a Monte 
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Carlo realization of a single stellar mass (1 .4M0) cluster with a flat distribution in stellar ages, and an H oc r/ scaling law for the likelihood of individual events. 



is illustrated in Figure 15, where, at low SMBH masses, the 
typical event is a MS disruption with evolved star disruptions 
contributing about 15%. As the black hole mass grows, an in- 
creasing fraction of MS stars are swallowed whole rather than 
disrupted, eventually leading to the sharp cutoff in the flaring 
rate seen in Figure 15. 

The limiting black hole mass scale of lO^M© coincides with 
a change in the SMBH flaring diet. In the right panel of Fig- 
ure 15, the typical disrupted star at high black hole masses is a 
giant branch star. At the very high mass end Mbh ^ lO^-^M©, 
and even HB stars are swallowed whole. The transition in 
the typical evolutionary state of disrupted stars as a function 
of black hole mass may also be inferred from Figure 16, in 
which the average tidal disruption mass feeding curve is plot- 
ted as a function of SMBH mass. For a population of SMBHs 
with Mbh ^ 1O^M0, the average M curve is dominated by MS 
disruptions with some perturbation at late times due to disrup- 
tions of post-MS stars. The representative M curves remain 
relatively unchanged until the black hole mass reaches the 



characteristic scale of lO^M©, where they transition rapidly 
to being dominated by post-MS disruptions. While post-MS 
M curves peak near the Eddington rate for low mass black 
holes, they fall well below the Eddington rate for the most 
massive SMBHs, as illustrated in the center panel of Figure 
16. Supplying mass near the Eddington rate of black holes 
> lO^M© via tidal disruption is not possible because of the 
transition to disruptions of stars exclusively in evolved stages. 

Our finding of a limiting SMBH mass scale is at odds with 
the results of Syer & Ulmer (1999) who compute rates of gi- 
ant consumption by SMBHs in Nuker galaxies and predict 
a high rate of disruption of giant stars in the most massive 
galaxies. This is because Syer & Ulmer (1999) consider a 
process in which the evolution of closely approaching stars 
results in mass transfer events with the black hole when they 
evolve off the MS. Contrary to their assumption, however, this 
process results in many weak encounters with the black hole 
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Fig. 16. — The left panel and center panels show the average of Monte Carlo realized flaring events in physical units and as a ratio to the Eddington rate, 
respectively. For simplicity, we take all events to be /3 = 1.5 encounters of a \AMq ZAMS mass star. We draw events along the star's lifetime using an 
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n oc r^' scaling assuming a flat distribution of stellar ages in the nuclear cluster, then scale to the nearest of our simulations in Figure 9. As black hole mass 
increases, the dominant flaring event goes from a MS disruption to a post-MS disruption, with the transition being particularly sharp around Mbh ~ 1O^M0. 
Post-MS disruptions feed gas to the black hole over longer durations than MS disruptions at correspondingly lower peak accretion rates. When normalized to 
the Eddington accretion rate of the black holes, Medd = O-O2eQ\(Mbh/lO^M0) Mqjy'^, this amounts to a precipitous drop in the luminosity of tidal disruption 
powered SMBH flaring as black hole mass increases. The right panel shows the duty cycle of flaring luminosity, taking the overall loss cone flux from Wang & 

Merritt (2004), N^6.5x 10'^ (Mbh/lO^Mo)"^/"^. Tidal disruption flares may be expected to contribute to the luminosity function of AGN harboring low mass 
black holes. However, they only contribute at very low luminosities relative to Eddington and with a low duty cycle in more massive SMBHs. 



(/3 <C 1) rather than a single disruptive encounter^. The sum 
of many encounters with the black hole may produce an inter- 
esting population of tidally-heated stars (Alexander & Morris 
2003; Alexander & Hopman 2003), but, is unlikely to produce 
luminous flares. 

5.2. Contribution of giant star disruptions to local 
low -luminosity AGN 

Having established a limiting SMBH mass scale for tidal 
disruption flaring events, we examine how this mass limit 
manifests itself in the cumulative, low-luminosity output of 
local AGN. Tidal disruption flaring events are observed to rise 
sharply and then decay over long timescales approximately 
following the widely discussed t~^^^ power law. The contri- 
bution of the tidal disruption of solar-type stars to the local 
AGN luminosity function was investigated by Milosavljevic 
et al. (2006) who found that the long decay tails of tidal dis- 
ruptions can contribute to the integrated luminosity of AGN 
in gas- starved galaxies. The exact value of the tidal dis- 
ruption rate as well as its evolution with redshift (which is 
rather uncertain; although see Merritt 2009) determines the 
time at which gas-mode feeding of active galaxies drops to 
sufficiently low levels for the stellar disruption component to 
take over (e.g. Heckman et al. 2004). For example, Milosavl- 
jevic et al. (2006) find that with a tidal disruption event of 

^ 6.5 X 10-4 (Mbh/lO6M0)"^^^ the associated flaring lu- 
minosity is likely to be responsible for at least 10% of local 
(z < 0.2) AGN activity. 

Observational AGN surveys have emphasized a picture of 
cosmic downsizing of black hole growth. The most massive 
SMBHs grew during the quasar era (Yu & Tremaine 2002), 
while the black holes believed to be actively growing now 
are much less massive, around lO^M© (Heckman et al. 2004). 
Our finding of a limiting mass scale for tidal disruption pow- 
ered flaring is consistent with the observation that most local, 
massive SMBHs are quiescent. Tantalizingly, the distribution 

^ The orbital period of trapped stars about the black hole, P ~ 10"^ years, 
is much less than even the post-MS stellar evolution timescale, Tg ~ 10^ 
years, resulting in very small per-orbit changes in stellar radius and impact 
parameter, /5. 



of local AGN-hosting galaxies is biased to about an order of 
magnitude lower SMBH mass than the general galaxy popu- 
lation (Greene & Ho 2007), with a decrease in activity around 
Mbh ~ lO^^M© that Hes coincident with the limiting SMBH 
mass scale for tidal disruption flares. 

The right panel of Figure 16 examines the duty cycle of 
tidal disruption powered flaring in local AGN. Despite their 
lower frequency of disruption, giant star flares make a sizable 
contribution to the average duty cycle at low luminosities be- 
cause their power law decay tails have higher normalization 
than those of MS stars (see Figure 9). The combined effects 
of the decrease in stellar consumption rate (Wang & Merritt 
2004), the precipitous drop in the flaring fraction (Figure 15) 
and the change in the demographics of disrupted stars (Fig- 
ures 15 and 16) all conspire to ensure that tidal disruption 
powered flares contribute well below Eddington and with a 
very low duty cycle (a small number in the active state) on the 
most massive SMBHs. 

The relative contributions of MS and evolved stars to the cu- 
mulative, tidal disruption-fueled accretion disk luminosity of 
local AGN depends on the total amount of stellar mass fed to 
SMBHs at rates below their corresponding Eddington limits. 
Super-Eddington fallback results in an inefficient conversion 
of fallback mass into light, since the accretion rate onto the 
black hole (and, correspondingly, the accretion disk luminos- 
ity) is capped at or near the Eddington limit^. Because their 
disruption results in highly super-Eddington feeding of lower 
mass SMBHs around lO^M©, MS disruptions are rather in- 
efficient at converting mass stripped from the disrupted stars 
into accretion luminosity. Giant stars disruptions, on the other 
hand, peak at lower accretion rates and feed more of their 
mass to the black hole below the Eddington limit. Figure 17 
shows the fractions of accretion disk luminosity contributed 
by stars at different evolutionary states. Giant star disruptions 
make up for in efficiency what they lose in disruption fre- 
quency when compared to MS stars and, thus, contribute sig- 

^ The super-Eddington fallback phase may, however, give rise to outflows 
(rather than accretion of mass with low radiative efficiency) and a conversion 
of mass to light at larger radius as the ejected debris expands (Strubbe & 
Quataert 2009, 2011; Kasen & Ramirez-Ruiz 2010) 
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Fig. 17. — The fractional distribution of total tidal disruption fed accre- 
tion disk light emitted, £^disk. calculated by a Monte Carlo realization of the 
mass fed to SMBHs below their Eddington limit, again for a IAMq star and 
/5 = 1.5 encounters. At lower SMBH masses, MS disruptions peak at highly 
super Eddington black hole feeding rates. As a result, MS disruptions inef- 
ficiently convert stellar mass loss into accretion luminosity. However, giant 
stars feed more of their mass to SMBHs below the Eddington limit, resulting 
in a more efficient conversion to disk luminosity. Thus, while post-MS dis- 
ruptions occupy about 15% of disruption events by rate, in terms of total disk 
emission, they contribute nearly 40% to the ensemble. As the SMBH mass in- 
creases to 1O^M0, both MS and post-MS disruptions are sub-Eddington, and 
the fractional contributions to fi'disk are similar to the event rate convolved 
with the mass loss per event (see Figures 15 and 6, respectively). 



nificantly to the overall power of tidal disruption-fueled AGN. 

5.3. Detection of giant star tidal disruption flares 

The primary characteristic that distinguishes the tidal dis- 
ruption flares of evolved stars from those of MS stars is their 
comparatively long timescales. Tidal disruption events are 
typically characterized in transient surveys by their late-time 
~ t~^/^ decay (see Gezari et al. 2009; van Velzen et al. 2011; 
Saxton et al. 2012, for some recent examples). However, as 
seen in Figure 9, We are unlikely to observe the characteristic 
power-law decay in the evolved cases due to their extremely 
long characteristic decay times. Our simulations show that the 
rise up to peak luminosity in giant star flares occurs over sim- 
ilar timescales as the decay in MS disruption flares. Figure 
1 8 illustrates that typical observing cadences will detect giant 
star disruptions during their rise to peak brightness. Shown in 
the post-MS category in Figure 1 8 are the average mass feed- 
ing rates of the RG I and HB disruption simulations. Tidal 
disruptions of stars at the tips of the giant branches (the RG 
II and AGB models) result in flares with such long timescales 
that they would be difficult to identify in transient surveys. 
The power law slope during the rise phase is quite steep, indi- 
cating that the brightness of a galactic nucleus would increase 
in a correlated fashion by many orders of magnitude over a 
period of a few years. 

AGN are infamous for exhibiting variability at a range of 
timescales, and one might worry that this would frustrate the 
detection of giant disruption flares. However, the power spec- 
trum of these variations is well below the amplitude expected 
from the disruption of a post-MS star. Van Velzen et al. (201 1) 
perform a detailed analysis of potential sources of tidal dis- 
ruption event confusion, including an empirical analysis of the 
amplitude of AGN flaring events in their Sloan Stripe 82 sam- 
ple of 1304 AGN. They find that the largest amplitude flaring 
events exhibit changes in flux of AF/F < 5, far less than the 
amplitude of the correlated rise expected from the disruption 
of an evolved star. More detailed studies of AGN variability 
also find small correlated increases (Kelly et al. 2009, 2011; 
MacLeod et al. 2010), giving credence to the idea that the tidal 
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Fig. 18. — The power-law slope inferred from two successive detections of 
hypothetical MS and post-MS tidal disruption flare transients given different 
observing cadences. Observing cadences that initially detect MS disruption 
events in the plateau or decay phases will catch the rise to peak in post-MS 
disruption flares. This plot assumes that flares above an intrinsic luminosity 
of 10^^ erg s~^ may be detected. This is roughly equivalent to the bolometric 
luminosity at which the recent tidal disruption flare PSl-lOjh at z = 0.1696 
was initially detected (Gezari et al. 2012). At the Eddington limit, we as- 
sume the luminosity of these flares plateaus but is convolved with a normally 
distributed variability factor with a = 0.01(Ar/lyr)^/^. The flares are all as- 
sumed to occur at an impact parameter of /3 = 1.5. The post-MS curve is the 
average of the RGI and HB models (see Figure 9). In both cases, the lines 
plotted are the median while the shading shows the ±lcr region for a Monte 
Carlo realization of detection times relative to the observing cadence. The 
inset plots show the detectable fraction of events, /det, defined as those typ- 
ically visible with at least two observations above the 10"^^ erg s~^ threshold 
given the cadence. 



disruption of an evolved star would imprint a unique signal in 
the AGN flaring history. 

The prospects for detecting giant star flares in the near fu- 
ture are, in fact, quite promising. While these events are less 
common than MS disruptions, they should be commonly de- 
tected in future transient surveys, like LSST, which are ex- 
pected to yield hundreds of MS disruption events per year 
(Strubbe & Quataert 2009). The insets in Figure 18 show that 
finding giant star disruption events does not require short ca- 
dence surveys. It is even possible that the signatures of giant 
star disruptions are lurking in current data sets and may be 
detectable through binning the data at lower cadences. Tidal 
disruption events with well sampled light curves that facilitate 
detailed comparison with theoretical M profiles offer the best 
hope of characterizing the physics of tidal disruption. 

5.4. Caveats and prospects 

In this work we have assumed that when the tidal radius is 
less than the Schwarzschild radius stars are swallowed whole 
without producing a flare. The impact parameter at which this 
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transition occurs, however, is likely sensitive to the orientation 
of the encounter relative to the spin of the black hole (Kesden 
2012; Haas et al. 2012). The exact fraction of events which 
are promptly swallowed or produce flares thus depends on a 
general relativistic description of the black hole's potential. 

A great deal of uncertainty lies in the properties of the nu- 
clear star clusters from which stars are fed into disruptive or- 
bits. We know very little about stellar kinematics in galac- 
tic centers other than our own (e.g. Genzel et al. 2010). Our 
work has made the usual assumptions of a spherical nuclear 
star cluster which feeds stars to the black hole by a two-body 
relaxation-driven random walk in angular momentum space. 
Theoretical work has emphasized a wide variety of effects 
which may enhance the rate at which stars are fed towards 
the black hole. A triaxial potential can feed stars to the loss 
cone collisionlessly through chaotic, boxy orbits (Magorrian 
& Tremaine 1999; Merritt & Poon 2004; Merritt & VasiHev 
2010). Rings or disks of stars, if present, would feed stars 
to the black hole at an enhanced rate through secular insta- 
bilities (Madigan et al. 2011). Another secular effect is reso- 
nant relaxation, which almost certainly enhances the loss cone 
flux to some extent for the most bound stars (Magorrian & 
Tremaine 1999; Hopman & Alexander 2006; Giirkan & Hop- 
man 2007). Finally, a second massive body, like an insprial- 
ing black hole or a giant molecular cloud could easily induce 
large angle scatterings of stars (Ivanov et al. 2005; Perets et al. 
2007; Chen et al. 2008, 2009, 2011; Wegg & Bode 2011). All 
of these processes result in more completely filled loss cone 
than the conservative scenario we've outlined here. 

Further, even the simplest characteristics of galactic center 
stellar populations, like the age and mass distributions of stars, 
are wildly uncertain. Stars are recycled from the galaxy in 
general as the interaction between tightly bound stars and the 
more weakly bound galactic stellar population lead to a time 
evolution of the nuclear cluster compactness and properties 
(Merritt 2009). In fact, our own galactic center is marked by 
a striking and enigmatic population of massive, young stars 
in the central ~ 0.1 pc (e.g. Schodel et al. 2007) that lies in 
contrast to the underlying, old stellar population. There are 
various alternative lines of evidence further suggesting that 
the initial mass function (IMF) of stars surrounding our own 
galactic center may be very different from that representative 
of the bulk of the star formation (Maness et al. 2007; Bartko 
et al. 2010). 

The dynamics of a cluster comprised of a spectrum of stel- 
lar masses is also expected to operate somewhat differently 
than that of a single mass cluster. In power law galactic cen- 
ters, which are presumed to be dynamically relaxed (Bahcall 
& Wolf 1976, 1977), one expects that massive stars, bina- 
ries, and stellar remnants to segregate to tightly bound orbits 
over the cluster's relaxation time. An illustrative example of 
these effects comes comes from simulations of globular clus- 
ters harboring intermediate mass black holes (IMBHs), which 
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have been performed with direct N-body integrations (e.g. 
Baumgardt et al. 2004a,b; Trenti et al. 2007) and Monte Carlo 
techniques (e.g. Giirkan et al. 2004). For example, Baumgardt 
et al. (2004b) consider a multi-mass stellar cluster and find 
that the massive stars segregate to form a tightly bound cusp 
around the IMBH, while the less massive stars relax to a shal- 
lower distribution. The scaling of the specific tidal disruption 
rate, /i, is modified by the varying radial profile from which 
stars are fed. In this case, more massive stars are fed into the 
loss cone preferentially as compared to low mass stars, and 
the rates of giant star disruption in their simulations are 15- 
20% of the total stellar disruption rate, a factor of ~ 2 greater 
than one derives from a simple weighting over the stellar mass 
spectrum. 

The detection of tidal disruptions holds great promise for 
revealing the properties of quiescent SMBHs and the nuclear 
star clusters that surround them, which, as we have discussed, 
are rather unconstrained. The universe is now being explored 
in a panchromatic way over a range of temporal scales, lead- 
ing toward a more complete and less biased understanding 
of its constituents. Flares from the disruption of stars in 
all evolutionary states will certainly be observed in the not- 
too-distant future. The recent observations of PSl-lOjh pre- 
sented by Gezari et al. (2012) offer a glimpse into disruptions 
of stars not on the MS. Taken in a statistical sense, the ob- 
served rates of tidal disruption and, in particular, the relative 
rates of disruptions of different stellar evolutionary stages, 
will hold tremendous distinguishing power in both the dy- 
namical mechanisms typically operating in galactic centers 
and the properties of the populations of stars themselves. Will 
the most common tidal disruption flares originate from young, 
massive stars in disks? Or from old stars fed from the near 
the edge of the black hole's sphere of influence? It is within 
the capacity of planned surveys to obtain lightcurves of gi- 
ant star tidal disruption flares that are well resolved tempo- 
rally. These lightcurves will challenge our understanding of 
the physics of tidal disruption, and through comparison with 
simulated events, offer a window into the extreme nature of 
close encounters between stars and SMBHs. 
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